Conceptual design of nanostructures for efficient photoinduced phase transitions 
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By means of Monte Carlo simulations on a kinetic model, we demonstrate that the efficiency of a 
photoinduced phase change can in general be enhanced drastically by using a superstructure of an 
appropriate combination of two components. This is due to the accelerated nucleation of converted 
domains in the structural blocks relatively close to local instability. The present mechanism provides 
a general guideline on the design of photocontrollable materials with potential applications for 
memory and storage devices. 
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Control of material phases by external stimuli, espe- 
cially by light, has been a subject of much interest in 
recent years in view of potential applications for memory 
and storage devices in the future. It requires a mate- 
rial exhibiting bistability between two distinct phases. 
In many cases the global bistability of a solid has a local 
origin while the photo-conversion between phases 

usually proceeds in a cooperative fashion . This phe- 
nomenon is therefore called a photoinduced phase transi- 
tion, in analogy with ordinary (thermally induced) phase 
transitions. 

The actual observation of the phenomenon has been 
limited within a small number of materials: one or a few 
examples from each class of materials, including spin- 
crossover complexes || , Prussian blue analogues , or- 
ganic charge-transfer complexes and conjugated poly- 
mers H . The scarcity of materials is presumably related 
to the requirement that the two phases must be energeti- 
cally close to each other, otherwise the metastable phase 
would return to the stable phase quickly. 

Here we propose a general scheme to overcome this dif- 
ficulty. Suppose we have a bistable material a in which 
the metastable state B is much higher in energy than the 
stable state A, so that photo-conversion is difficult. Let 
us then prepare a companion material [3 in which the 
relative stability is opposite (see Fig. 1(a)). If we com- 
bine the two materials in an appropriate fashion, e.g., 
in a superlattice as shown in Fig. 1(b), and if there are 
sufficient cooperative interactions that favor the a and (3 
units being in the same state (either A or B), then the 
energy of the A phase in the total system can be made 
degenerate with that of the B phase. This is rather an ob- 
vious way of designing new suitable materials in general. 
However, this approach yields more than that: We will 
demonstrate in the following that the photo-conversion 
efficiently is drastically enhanced in such mixed struc- 
tures compared with uniform structures. 

We consider a phcnomcnological model described by 
the Ising Hamiltonian 0,|8| 



H = - 



(1) 



(i,j)=n.n. 



on a three-dimensional lattice, where the 'spin' variable 



Si = ±1 denotes the two states (B and A) of the site 
i with the energy difference 2si. The nearest neighbor 
(n.n.) coupling Jij is assumed to be a positive constant 
J, so that the neighboring sites prefer to be in the same 
state. We use a simple cubic lattice of N x N x N sites 
with periodic boundary conditions. N — 60 is a typical 
size used in our simulations. We consider the superstruc- 
tures consisting of two kinds of blocks, a and /?, alter- 
nating in one-, two-, or three-dimensional (ID, 2D, 3D) 
directions as shown in Fig. |l](b). The size of each block 
is n in the confined directions. (The system size N must 
be an integer multiple of 2n.) The site potential Si is 
assumed to be constant in each of the blocks, —A/2 and 
A/2 in the a and (3 blocks, respectively, corresponding 
to the situation illustrated in Fig. |l|(a). Note that the 
system is symmetric with respect to the interchange of A 
and B (a and (3). The case of A = corresponds to the 
uniform system. 

For A > 0, the a and (3 blocks favor the A and B 
states, respectively. However, the system still tends to 
be in a uniform state without phase separation, if the 
block thickness n is not very large and the mier-block 
interface energy proportional to J outweighs the intra- 
block energy gain proportional to |A|. The condition for 
uniformity is written as |A| < A c = AdJ/n, where d is 
the dimensionality of the superstructure. Under this con- 
dition, the uniform ground states are doubly degenerate: 
The phase of all the sites being in the A state has the 
same energy as that of all the sites in the B state. 

We have studied the kinetics of the model at finite 
temperatures under external excitations by means of a 
classical Monte Carlo simulation ||, wherein the spin 
flipping rate at each site j is the sum of the thermal and 
external terms: uij = wj + wj . When we consider selec- 
tive excitation from A to B, we set wJ(A — > B) = W and 
wj(B -»• A) = 0. In each Monte Carlo step (MCS), both 
thermal and external spin-flip trials are performed at ev- 
ery site once on average. The heat bath model is used for 
the thermal flipping rate wj , although we have checked 
that other algorithms such as the Glauber, Metropolis, 
and Arrhenius models give essentially the same results. 

Suppose the system is initially in the A phase and de- 
velops under an excitation rate W from A to B at tem- 
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perature k^T — J. Figure g(a) shows the fraction of 
sites in the B state as a function of the excitation time 
t for various W in the homogeneous case of A = 0. For 
weak excitation with W — 0.09, the system remains per- 
manently in the A phase, whereas conversion to the B 
phase occurs for larger W. The converted fraction does 
not increase linearly with time but grows rather suddenly 
after a certain time. This is a clear reflection of the coop- 
erative nature of the phase change. This type of behavior 
was actually observed in a spin-crossover complex, where 
the transition occurs with the absorbed photon flux of 
about 10 18 cm^s" 1 §. 

We move on to the case of a superstructure with A ^ 0. 
Figure ||(b) displays the time evolution for a weak excita- 
tion rate W = 0.05 at k^T = J in the 2D superstructure 
with the block size n = 1 and A = 0.8A C = 6.4J. The 
phase change occurs in a short time in contrast to the 
case of A = 0, for which there is no indication of a phase 
transition. (Note the different time scales of Fig. 2(a) and 
Fig.2(b).) 

To clarify the underlying mechanism of the accelerated 
phase transition, we plot in Fig. ^(b) the converted frac- 
tions in the a blocks and in the f3 blocks separately. The 
converted fraction of the (3 blocks grows quickly from the 
beginning, and then that of the a blocks gradually fol- 
lows. By tracking the time evolution of all the sites, it 
has been found that a nucleated domain of B states in a 
(3 block does not decay so quickly as in the case of the 
uniform structure because of the local stabilization en- 
ergy A, thereby allowing further growth of the B domain 
around the nucleus firstly in the f3 block and then onto 
the surrounding a blocks. Although the energy differ- 
ence A counteracts as an energy cost in the a blocks, 
the growth process continues once the stable nucleus is 
formed. 

The time of continuous excitation required for induc- 
ing the phase transition, t p , depends on the excitation 
intensity W. The IF-dependence of t p is shown in Fig. || 
for various A in the 2D superstructure. Here t p is de- 
fined by the time of the converted fraction exceeding 0.8. 
There is a threshold W c below which the phase transi- 
tion does not occur for each A, and t p sharply diverges 
from above W c . We have confirmed the critical behavior 
by extending simulations up to 100,000 MCSs. Similar 
results have been obtained for ID and 3D superstruc- 
tures. The existence of a threshold light intensity for the 
photoinducedphase transition was also suggested in real 



at which t p diverges 



iced p 

materials |Q. The threshold W t 
is reduced drastically with increasing |A|. The results in 
Fig. H imply that the required time t p at each W dras- 
tically decreases with increasing |A|. For example, at 
W—0.12, the phase transition in the case of A = 0.9A C 
requires only one tenth of the excitation time for the ho- 
mogeneous structure with the same excitation intensity. 

In conclusion, by using nanostructures, we can dras- 
tically reduce the threshold light intensity and enhance 
the efficiency of phase switching in both the forward and 
backward directions. With this principle, which does not 



rely on the details of materials, the range of useful materi- 
als can be extended substantially, because an appropriate 
combination of inferior materials can produce a superior 
material for photoinduced switching. 

The present mechanism was hinted by the study of 
the photoinduced magnetism in cobalt-iron cyanides 
]l|,pyTl]-[T^] , which can be considered as a mixed system 
of different local environments with cooperative interac- 
tions due to lattice distortions Q . It would be very in- 
teresting to make a superstructure in this class of materi- 
als as well as in other materials such as perovskite man- 
ganese oxides |l^|l6|], for which a photoinduced metal- 
insulator transition was observed in the bulk |T?j ] . 
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FIG. 1. (a) Schematic bistable potentials with opposite 
relative stabilities in a and j3 units, (b) Schematic one-, two-, 
three-dimensional superlattice structures consisting of the a 
and B blocks. 
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FIG. 2. Monte Carlo simulation on the kinetics of the ex- 
citation-induced phase transition in the system of 60 x 60 x 60 
sites at feT = J- (a) The fraction of converted sites as 
a function of the excitation time t in Monte Carlo steps 
(MCS) with A = for various excitation rate W per MCS. 
(b) The converted fraction in the 2D superstructure with n 
= 1 for W = 0.05 per MCS at k B T = J, in the case of 
A = 0.8A C = 6.4J (thick line) compared with the case of 
A = (thin line). Solid, broken and dotted thick lines show 
the total fraction and the partial fractions of the a and f3 
blocks, respectively. 
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FIG. 3. Dependence of the conversion time t p on the 
excitation intensity W in the 2D superstructure with n = 1 
at ksT = J for various A with A c = 8 J. 
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